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Abstract 



The connection between scarring and tunneling in chaotic double-well po- 
tentials is studied in detail through the distribution of level splittings. The 
mean level splitting is found to have oscillations as a function of energy, as 
expected if scarring plays a role in determining the size of the splittings, and 
the spacing between peaks is observed to be periodic of period 2ttTI in action. 
Moreover, the size of the oscillations is directly correlated with the strength 
of scarring. These results are interpreted within the theoretical framework 
of Creagh and Whelan. The semiclassical limit and finite-?i effects are dis- 
cussed, and connections are made with reaction rates and resonance widths 
in metastable wells. 
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I. INTRODUCTION 



Many chemical reactions must proceed through a potential barrier before the final dis- 
sociation of the products of the reaction can take place. Thus, the reaction rate is governed 
by tunneling |]. Radioactive decay of nuclei also involves crossing a potential barrier, dic- 
tated by a combination of short-ranged strong binding forces and a longer-ranged Coulomb 
repulsion 0. A typical experimental situation is the absorption of a slow neutron; the re- 
sulting metastable nucleus decays predominantly by tunneling through a saddle point, and 
the distribution of resonance widths is well described by random matrix theory ||. Another 
important example of quantum tunneling is the conductance of mesoscopic devices [f|]. In 
more than one dimension, when the classical dynamics can be chaotic, quantum chaos plays 
a role. For instance, in the tunneling diode junction, in which electrons are driven by an 
applied electric field but must tunnel through potential barriers at either end of the device, 
the dynamics, in the presence of a magnetic field, is chaotic and scarring of a short periodic 
orbit is known to dominate the conduction [Q. In this paper we study chaotic tunneling 
in a simple model potential, and establish the connection to scar theory, both qualitatively 
and quantitatively. We prefer to study tunneling by calculating splittings in a double-well 
potential rather than resonance widths in a metastable well, because the former is an easier 
computational task, as we will see later in Sec. [TT]. The relation between splittings in a 
double-well potential and resonance widths in a metastable well will be discussed imme- 
diately below, where we will observe that the two quantities are related by known overall 
normalization factors. 

The double-well potential we will work with, in two dimensions, has the very simple form 
V(x, y) = x*- x 2 x 2 + ay 2 + Xx 2 y 2 + ]T b^^+^-y,?)!^ . (1) 

The parameters xo, a, A, fej, Xi, yi and <Tj will be specified below. For x\ > 0, a barrier along 
the |/-axis separates the potential at low energies (E < 0) into two wells. The Xx 2 y 2 term 
ensures that the potential is not separable, while the Gaussian perturbations, the parameters 
of which may be changed at will, allow us to generate an ensemble of statistically independent 
eigenstates near any given energy. If the positions and heights of the Gaussians are chosen 
suitably, the potential has reflection symmetries in the x and y directions. The symmetry 
under reflection in y means that there will always be a short periodic orbit on the x-axis, in 
each of the two wells. 

The principal result of this paper is that the size of the level splittings in the two- 
dimensional double well, the classical dynamics of which is chaotic, is directly correlated 
with the scarring of eigenfunctions along the x-axis, which we believe to be the primary 
channel for tunneling. The evidence for this is that (1) the distribution of splittings, once 
the average exponential trend in energy is scaled out, displays oscillations as a function 
of energy of the sort expected by scar theory. In particular, the action distance between 
successive peaks is precisely h. However, the dependence of the mean splitting on action is 
not well reproduced quantitatively by linear scar theory, because non-linear effects appear 
to be important. (2) The rescaled splittings are correlated with the overlap of the eigenstate 
with a Gaussian test state lying on the periodic orbit on the x-axis, a measure of scarring. 
(3) The correlation of the splittings with the overlaps is further supported by the fact that 
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the distribution of overlaps has the same energy- dependent oscillations as the distribution 
of splittings. 

According to one-dimensional WKB theory, the splitting AE = -E a nti-s y mm — -Esymm in a 
symmetric double well is given in the semiclassical limit by 



(2) 




while the resonance width T for the state at the same energy in a metastable well is 

'M (3) 

Here ui is the frequency of the classical periodic motion at energy E symm and the imaginary 
action for going under the barrier is 

S= / dx Jv(x) — E , (4) 

•J—Xtp 

x tp > being the position of the classical turning point at energy E syram ||. In order for the 
semiclassical theory to apply, E symm must be sufficiently far below the barrier that S/Ti~^> 1. 
Note the factor of two in the exponent in Eq. (|3]) . 

These one-dimensional formulae may be generalized to the two-dimensional (or higher- 
dimensional) potential well as follows. We expect tunneling to be dominated by paths that 
cross the barrier close to the x-axis, which has the smallest action integral S. The splitting 
AE will then be proportional to the exponential factor e~ s / h . 

The correct generalization of the frequency of attempts to cross the barrier, d, to two 
or more dimensions, is the frequency with which one returns to a Planck-sized cell in phase 
space that lies on the horizontal periodic orbit. This horizontal periodic orbit is the real 
continuation of the least-action path across the barrier. The time for returning to such a cell 
(or to any other cell in an ergodic well) is the Heisenberg time Tjj = h/A(E), where A(E) is 
the mean level spacing near energy E (i.e., the spacing between doublets in the double- well 
system). Then the frequency of attempts to cross the barrier is just proportional to the mean 
level spacing A(E). Thus, we expect on general grounds that the splitting AE is given in 
order of mae nitude by AE ~ A(E)e- s( - E)/h , which gives us the trend of the splittings 
function of energy. This expression for the mean splitting will be confirmed by the exact 
semiclassical theory to be discussed below. 

For any given state, we should expect that its splitting will be large or small compared 
with the mean value at that energy according to whether its amplitude is large or small on 
the horizontal periodic orbit which leads to optimal tunneling. For simplicity we can study 
the wave function amplitude near the turning point of the horizontal periodic orbit. The 
value of the wave function at the turning point in the two-dimensional chaotic system is 
(ignoring scar-related effects which are the main focus of this paper) given approximately 
by a Gaussian-distributed random variable, as random matrix theory would predict. Thus, 
|?H^tp)| 2 has, according to random matrix theory, a Porter-Thomas distribution for all 
energies far enough below the barrier (near zero energy the horizontal periodic orbit becomes 
stable and the distribution of splittings rolls over to one having many more large and small 
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splittings, corresponding to wave functions that live near or avoid this stable orbit). However, 
what is relevant to tunneling in d dimensions is not just the value of the wave function exactly 
at the turning point but rather its behavior in a whole ft d_1 — sized region surrounding the 
periodic orbit. We shall see below that the right quantity to consider is the inner product 
of the wave function with a Gaussian centered on the periodic orbit (at the turning point or 
at some other location). This has as well a Porter-Thomas distribution, within the random 
matrix theory approximation. Scar-related effects and finite- ft effects on the distribution of 
splittings will be discussed below. 

The viewpoint summarized above is in agreement with the theoretical work of Creagh 
and Whelan [[/J. First, they find that the mean splitting (AE) at a given energy is given 
by the product of an exponential factor e~ s ^ n , a factor proportional to the mean spacing 
between doublets, as described above, and a third factor that carries information about the 
monodromy matrix of the (imaginary time) tunneling orbit. Then they show that, for chaotic 
and symmetric double wells, the splitting for a particular eigenvalue, relative to the mean 
splitting, may be written in the semiclassical limit as a matrix element of the wave function ip 
near the real continuation 1Z of the complex trajectory that passes through the barrier with 
minimum (imaginary) action. This imaginary-time trajectory and its real-time continuation 
may be thought of as the optimal route for tunneling. The matrix element AE ~ (ip\T\ip) 
involves integration over a Poincare surface transverse to the real continuation 1Z. The 
kernel T is a semiclassical Green's function which, in the approximation that the dominant 
contribution to the tunneling matrix element comes from the neighborhood of 7Z, becomes 
a Gaussian centered on the intersection of 1Z with the Poincare section. The width of the 
Gaussian is of 0(ft 1//2 ) in both directions (e.g. y and p y ) tangent to the surface of section. 
[The results may of course be easily generalized to dimensions d > 2, where the resulting 
Gaussian has width of 0(ft 1//2 ) in all 2d — 2 directions along the surface of section.] 

In the case that the real continuation 1Z happens to lie on a short periodic orbit, which will 
always be true when a reflection symmetry across the x-axis is present, this matrix element 
may be regarded as a measure of scarring on the periodic orbit. The Creagh- Whelan theory 
predicts, therefore, that strong scarring should be correlated with large splittings, confirming 
the intuitive expectation that high tunneling rates should occur for those wave functions 
that have large amplitude along the path with optimal tunneling. Neither in Ref. J7| nor in 
Ref. P|, however, do the authors demonstrate the connection between scarring and splittings 
on a state-by-state basis. In the latter work, they confirm their formula for the tunneling 
matrix element by deriving from it an analytical prediction for the statistical distribution 
of splittings. This prediction is in agreement with numerical calculations for potentials in 
which the real continuation of the optimal tunneling orbit is not a periodic orbit; when 
it is, the random-matrix assumption in their derivation breaks down due to scarring on 
this periodic orbit. Thus, the present paper, while confirming the predictions of Creagh and 
Whelan, goes beyond their results by establishing conclusively the link between scarring and 
tunneling and by showing, with better statistics, that the distribution of scaled splittings 
indeed becomes approximately Porter-Thomas, but only after scarring effects have been 
taken into account. 
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II. METHOD 



The wave functions and splittings were calculated numerically using the discrete variable 
representation ||. The matrix elements of the position operators X and Y and of the kinetic 
energy operators K x and K y were first evaluated analytically using standard identities, in 
a basis of up to the first 300 Gauss-Hermite functions in each dimension. Then, in order to 
take advantage of the two reflection symmetries in x and y, the two operators X 2 and Y 2 
were diagonalized. The reason for using X 2 and Y 2 instead of the usual choice of X and Y 
is that X 2 , K x , Y 2 and K y are all block-diagonal, connecting only basis elements within one 
of the four symmetry classes (even-even, even-odd, odd-even and odd-odd). Since we are 
interested only in even-even potentials of the form V(X, Y) = YU fi{X 2 )gi{Y 2 ) we can just as 
well compute V at the eigenvalues of X 2 and Y 2 as at those of X and Y, but using the basis 
obtained by diagonalizing X 2 and Y 2 ensures that the final Hamiltonian H = K x + K y + V 
will be itself also block diagonal. The four symmetry classes may therefore be analyzed 
separately. In the basis chosen the potential V is of course diagonal, while, in two or more 
dimensions, the kinetic energy matrix will be sparse. More precisely, if N is the dimension 
of one of the blocks in the Hamiltonian matrix, the total number of non-zero matrix entries 
scales as N 1+1 ^ d , or N 3 ^ 2 in the two-dimensional system. Since we require large values of N 
in order to observe semiclassical behavior, a sparse matrix routine is the method of choice 
for diagonalizing H . The accuracy of the computed eigenvalues was tested for convergence 
under increase of N, and for the parameters given below we found convergence to ±10~ 12 
for N fa 3500, corresponding to about 200 Gauss-Hermite functions in the x direction and 
about 100 in the y direction. 

The amount of phase space covered by the region E < 0, and thus the number of states 
under the barrier and the computation time, increases very rapidly with xq. If all other 
parameters in the Hamiltonian are kept fixed, the number of states grows as X®, and so the 
largest value of Xo we can easily attain is about xq = 6, for a = 1 and A = 10. At these 
parameter values, each well has a depth of Xq/4 = 324 and about 100 bound states, where 
h is taken to be unity here and in the following. The typical level spacing near the top 
of the well is ~ 1, and the splittings range from < 1 near the top of the well to < 10~ 6 
near E = —50; below this energy many of the splittings become too small (< 10~ 12 ) to be 
resolved numerically. Therefore we take E = —50 to be the lower cutoff for the energies to 
be analyzed in Sec. |Tj. 

III. RESULTS 

The potential given by Eq. (|1|), apart from the Gaussian perturbation, is mostly in- 
tegrable for all energies except those near the top of the barrier. When the Gaussian 
perturbations are introduced, the classical mechanics becomes more chaotic, but if these 
perturbations are too small it is still possible that they would not be seen by the quantum 
mechanics, which would remain effectively integrable. Thus, in order to render the quantum 
mechanics at energies corresponding to the bound states chaotic as well, it is necessary to 
introduce a Gaussian perturbation that is at least as large as the wavelength in question, and 
whose height is comparable in magnitude to the depth of the potential well. The simplest 
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choice is to place large Gaussian perturbations above the minima of the potential well at 
(±xo/v2, 0), which will be seen by every bound trajectory as it crosses the center of the 
well and which thus effectively make the dynamics chaotic at energies down to the lowest 
considered (E = —50). This was checked classically by examining the Poincare surfaces 
of section and may also be seen to be true quantum-mechanically in Fig. [I], where typical 
eigenfunctions are shown. Here, we have chosen the parameters of the double-well potential 
to be Xq = 6, a = 1 and A = 10, and for the central Gaussian perturbation we use a height 
b = 150 (to be compared with a well depth of V(±xq/\/2, 0) = —324) and a width a = 0.5. 

We generate an ensemble of 625 systems by placing four further Gaussians (and their 
reflections in x and y) at x — ±2, ±3, ±4, ±5 and y = ±1, with heights hi = 20rii, 
ni,...,7i4 = 1,...,5, and with equal widths Oi = 0.5 as above. A contour plot of the 
potential for a typical member of the ensemble is given in Fig. 0; note the symmetrical 
distortion of the contours due to the perturbation. 

We proceed to analyze statistically the splittings between states in the even-even and 
even-odd sectors. The results, for the parameters described, are given in Figs. ^| and |j. As 
expected, the size of the splittings falls off exponentially with decreasing energy in Fig. |3], as 
the barrier becomes wider and tunneling is suppressed. The trend is approximately linear 
on a semi-log plot, over six decades as E varies from to —50. In Fig. |4] we rescale the 
splittings as a function of energy by s — > s/e~ s( - E \ We find very pronounced oscillations 
in the distribution of splittings as a function of energy. As discussed in Sec. |, we expect 
theoretically that the rescaled splitting should be proportional to the overlap in a Poincare 
section of the eigenfunction with a Gaussian on the horizontal periodic orbit. If the Gaussian 
may be assumed to have area exactly h in the Poincare section, such overlaps are described by 
scar theory |lTi|-|T2| . We expect the results to be qualitatively the same even if the Gaussian 



in the Poincare section given by the theory of Creagh and Whelan is not a minimum- 
uncertainty state, as long as it is not too large compared to h (in the latter limit, scar effects 
must go to zero in accord with the Schnirelman ergodicity theorem [TB|]). In fact, in the 
data presented below the area of the Creagh-Whelan Gaussian ranged from 1.5h to 4h. The 
prediction of scar theory is that, at a given energy, the distribution of splittings should be 
Porter-Thomas, and that the mean wave function intensity and therefore the mean splitting 
should oscillate as energy is varied by an amount that depends on the Lyapunov exponent 
and the monodromy matrix of the unstable periodic orbit. An important confirmation of 
the scarring picture is obtained when we plot, in Fig. [5j, the rescaled splittings versus the 
action (divided by 2ir) of the horizontal periodic orbit at the energy eigenvalue. We find 
that the oscillations are periodic in action with period 2tt, which indicates that the EBK 
quantization condition for scarring holds. The EBK quantization condition for the action 
reads A = 2n(n + l/2+n c /4) where n is an integer and n c is the number of conjugate points 
in one period of the orbit (n c = 3 in the case of the horizontal orbit in our system). 

A direct correlation between splittings and scarring is found by plotting, in Fig. |6|, the 
rescaled splitting of each eigenvalue versus the overlap of the corresponding eigenfunction 
with a Gaussian test state lying on the horizontal periodic orbit, a measure of the degree 
with which this eigenfunction is scarred. The two quantities are correlated, with a slope of 
2 on a log-log scale. The correlation coefficient of the logarithms is 0.78; it may be that 
the degree of correlation would be improved if the Gaussian were chosen to be properly 
aligned with respect to the monodromy matrix of the optimal tunneling path. The observed 



6 



correlation nevertheless confirms that there is a direct connection, on a state-by-state basis, 
between scarring and tunneling, as predicted by the theory of Creagh and Whelan 0. As a 
check on our results, we show in Fig. that the overlaps display the same energy- dependent 
oscillations as do the splittings, as they must if the phenomenon of scarring underlies the 
behavior of both. 

The connection between scarring and tunneling can be tested quantitatively in two ways. 
First, scar theory in the semiclassical limit predicts that the short-time (smooth) envelope 
describing the oscillations in the mean rescaled splitting versus action is given by the Fourier 
transform of the autocorrelation function A(m) = (<j)\(f)(m)) , where is a Gaussian wave 
packet (living in the Poincare section) centered on the horizontal periodic orbit and <j)(m) 
is its iterate after m bounces. The Gaussian (f) is chosen to have the same orientation and 
aspect ratio in the (y,p y )-plane as the Gaussian called for by the Creagh- Whelan theory, 
but linearly rescaled so as to have area h as needed for scar theory. Linearizing the dynamics 
around the horizontal periodic orbit we find, when the Gaussian wave packet is optimally 
aligned along the stable and unstable manifolds of the periodic orbit, 

A(m) = 1 , (5) 
V cosh Am 

where A is the Lyapunov exponent of the periodic orbit. In the special case of orthogonal sta- 
ble and unstable manifolds, a circular Gaussian will be one example of an optimally aligned 
wave packet. If the Gaussian is not optimally aligned, this formula may be generalized as 
follows: 



A(m) 



det(M) 



\ det(M+ (J~ m ) T MJ 



where M describes a Gaussian of the form (const) exp(— x T Mx) with x = (y,p y ) T represent- 
ing the coordinates in the surface of section, and J is the Jacobian of the Poincare mapping 
evaluated at the periodic orbit (Tr J = 2 cosh A). The matrix M is computed as specified in 
the Creagh- Whelan theory from the monodromy matrix of the complex orbit that begins at 
the Poincare section on the right, goes through the barrier, and ends at the Poincare section 
in the left well J7J . Here, the Lyapunov exponent A and Jacobian J vary with energy over the 
range —50 < E < —9. At higher energies, the trajectory spends less time near the Gaussian 
at the center of the well, and thus experiences less deflection, leading to greater stability, 
eventually becoming stable for E > —9. The short-time envelope obtained as the Fourier 
transform of A(m) in either Eq. (^) or Eq. (Q) may be compared with the mean rescaled 
splitting plotted versus action. As shown in Fig. with either form of the autocorrelation 
function we do find peaks in the predicted envelope of splittings at the right values of action 
for energies E < —9 (for energies E > —9 the horizontal periodic orbit becomes stable, so 
the scar theory does not apply and no prediction about the distribution of splittings can be 
made), but the heights of the maxima and minima between the peaks are not well repro- 
duced. The contrast predicted by Eq. (JH) is closer to the numerical data than that predicted 
by Eq. (|5]). The quantitative failure of semiclassical scar theory is attributable to the fact 
that, for our parameter values, the linearizable region around the horizontal periodic orbit 
is not large compared to h. In fact, the size of the linearizable region is only about 0.15/i 
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for the energies considered. Its size, however, is approximately independent of energy for 
—50 < E < 0, and this may explain the weak dependence of the peak height on energy 
observed in the numerical data. Nevertheless, we see that not only is scarring associated 
with larger splittings in the coarse sense of Fig. ^ but also the enhancement factor in the 
distribution of splittings, as a function of energy or of action, does oscillate with energy or 
action in agreement with the analytical prediction of scar theory; only the precise magnitude 
of these oscillations remains unexplained within the present linear theory. 

A second quantitative test of scar theory in relation to tunneling is to examine the change 
in the distribution of splittings upon change in the Lyapunov exponent. The horizontal 
periodic orbit can easily be made more stable by keeping the height of the main Gaussian 
bump fixed at 150 while increasing its width. An ensemble of eigenstates and associated 
splittings was computed, just as above, for a larger value of the bump width, namely 0.63 
instead of 0.50. The expectation from scar theory would be for the distribution of splittings 
to have many more smaller and larger splittings at the resulting smaller Lyapunov exponent. 
At a = 0.5, the horizontal periodic orbit is stable down to E — —9.1; the Lyapunov exponent 
then increases from zero with decreasing energy to a value of A = 2.0 at E = —50. For 
a = 0.63 it is stable down to E = —49.8 and attains only a value of A = 0.11 at E = —50. 
The numerical data, however, show no marked difference between the two computations at 
a = 0.5 and o = 0.63; see Fig. [| The lack of any difference between the distributions of 
splittings despite the difference in stability is an indication that we are not far enough into the 
semiclassical limit (see discussion below). In both cases the distribution of rescaled splittings 
(see the histograms in Fig. |9|) has many more small and large splittings, and consequently 
fewer splittings around s/(s) = 1, than a Porter-Thomas distribution would have (except 
for the sharp cutoff at s/ (s) > 5, which will be discussed below). Thus, the prediction of scar 
theory that there should be many more small and large splittings, relative to the prediction 
of random matrix theory, is confirmed. Also, the divergence of the probability distribution 
near zero splitting in the case of scarring on the real continuation of the optimal tunneling 
path differs markedly from the results, both analytical and numerical, of Creagh and Whelan 
|| for the case when the real continuation is not a periodic orbit, which show a probability 
distribution tending to zero at zero splitting. Our numerical results for the scarring case 
improve on their statistics and allow us to discern the scar corrections to Porter-Thomas 
behavior. In particular, we note the excess of very small splittings; these correspond to the 
phenomenon of antiscarring, as seen in Fig. || at actions halfway between values of action 



given by the EBK quantization condition for maximal scarring. As studied by Kaplan WA 



in an open quantum system coupled to the environment by one channel located on a short 
unstable periodic orbit, antiscarring causes the probability to remain in the system at times 
large compared to the Heisenberg time to be substantially enhanced relative to the prediction 
of random matrix theory. Therefore, we must expect that antiscarring, which we have 
demonstrated now for the case of level splittings in a smooth chaotic double-well potential, 
would markedly alter, away from random matrix theory predictions, the distribution of 
resonance widths in a chaotic metastable potential, and also the long-time probability to 
remain in such a well. 

In Fig. [10] we show the relation between the rescaled splitting AE/e~ s and the mean 
level spacing A, which decreases from about 3 at E = —50 to about 1 at E = 0, aside from 
some fluctuations. There is no direct correlation between AE/e~ s and A, thus refuting 
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the intuitive expectation that the tunneling rate should be proportional only to the rate of 
attempts to cross the barrier given by the classical motion, as discussed above in Section |. 
In the presence of scarring on the horizontal periodic orbit, tunneling is enhanced by the 
tendency to remain near the horizontal periodic orbit. At energies for which scarring takes 
place, the typical wave function intensity measured using a Gaussian at the turning point 
will be enhanced by a factor of 0(1/ A) compared with the naive expectation, where A is the 
Lyapunov exponent. At energies for which antiscarring, a tendency to avoid the horizontal 
periodic orbit, takes place, this typical intensity will be strongly suppressed, by an amount 
that is exponentially small in A for small A. The actual distribution of the rescaled splitting 
AE/e~ s versus the level spacing A(E) includes the same energy-dependent oscillations seen 
in Fig. as a function of A(E) rather than of E itself. It is evident, then, that chaotic 
tunneling in two dimensions must be thought of as a quantum- coherent phenomenon, in 
which the probability of tunneling through the barrier is greater if one comes back in phase 
when making repeated attempts to cross the barrier, as happens for scarred eigenfunctions. 
We also note that the horizontal periodic orbit becomes more unstable at lower energies, 
leading to smaller scar peaks in the mean wave function intensity on the orbit, and thus 
compensating to some extent for an increase in the mean level spacing at lower energies. 
This may partly explain the absence a clear trend in the data of Fig. [It]. 

We now discuss how our data are limited by the fact that we must work at finite h. 
First, the sharp cutoff at large splittings in the numerical data relative to the Porter- 
Thomas distribution in Fig. |9| is a finite-^ effect. This can be understood as follows. Let the 
Poincare surface of section have area N in units of h. Then the expected squared overlap 
(s) of an eigenstate with a Gaussian test state will be l/N, because the test state covers an 
area h in phase space while the eigenstate is, on average, spread evenly over the entire phase 
space. Now the cutoff arises from the fact that no matter how scarred or otherwise localized 
the eigenstate is, its overlap with a test state cannot be greater than unity. So s < 1 by 
construction, or s/(s) < N. Thus the cutoff increases to infinity in the semiclassical limit 
(i.e., as h tends to zero), even while (s) itself is decreasing. Assuming random matrix theory, 
the modified form of the Porter-Thomas distribution for finite N can be computed. One 
takes an ensemble of randomly oriented vectors in N dimensions, normalizes them so they 
lie on the unit sphere, and takes N times the square of the z-component. This quantity 
has mean 1 and a sharp cutoff at N. The Porter-Thomas distribution is recovered in the 
limit N — > oo. For an analytical form for the Porter-Thomas distribution for finite N, see 
Brody et al. in Ref. [|T5[|, especially their Eq. (7.10). In Fig. |TT] we see that the modified 
Porter-Thomas distribution for N = 6, corresponding roughly to the effective dimension 
of our Hilbert space, reproduces the cutoff in the numerical data of Fig. The scarring 
corrections (extra splittings at large and small s/ (s) with fewer splittings around s/(s) = 1) 



relative to the modified Porter-Thomas distribution are still present in Fig. [TT. 

A second test of the effect of finite H is to repeat the calculation at different values of 
h. Since we are near the computational limit already, we consider only the case of larger h. 
This is done by scaling the coordinates (x,y) — > (x',y') = (cx,cy), < c < 1. Under this 
transformation the potential becomes 



/-»rt /"'■^ f%& (~*^ 
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while the kinetic energy remains 

(8) 

dx' 2 dy' 2 y ' 

since the momenta are not affected by the transformation. The complete transformation of 
the Hamiltonian may be regarded as the product of three transformations: (i) the scaling 
of coordinates by a factor of c 1//2 and momenta by a factor of c _1//2 , which does not change 
the quantum mechanics, (ii) scaling both coordinates and momenta by a common factor of 
c 1 / 2 while also replacing the Hamiltonian H by cH, which preserves the classical mechanics 
exactly but is not area-preserving, and thus affects the quantum mechanics by changing 
the effective value of h, and (iii) scaling the Hamiltonian by a factor of 1/c, which trivially 
rescales the spectrum back into the original range. The reason we use this transformation 
is to keep the classical mechanics, all the periodic orbits, their stability properties, etc. 
unchanged as we change the effective value of h, so the results for different values of the 
effective H (which scales as 1/c) are directly comparable. 

For c = 0.8 we find the same oscillations observed previously in the distribution of 
rescaled splittings as a function of energy, only now there are four peaks in the range from 
E = —50 to E = compared to the five that we saw before, corresponding to a larger 
effective value of H in the new system. The distribution of splittings is given by the scarring 
corrections to the modified Porter-Thomas distribution for N = A now, compared to N = 6 
above. Thus, the same conclusions continue to hold but with the expected modifications for 
larger h. This indicates that at c = 1 we are far enough into the semiclassical regime to see 
characteristic semiclassical behavior for the locations of the scarring peaks, if not for their 
precise heights. 

IV. CONCLUSIONS 

We have demonstrated that scarring on the real continuation of the optimal tunneling 
path, if it is an unstable periodic orbit, enhances tunneling and thus leads to larger splittings 
between the symmetric and antisymmetric in x eigenfunctions at energies near the scarring 
energies (likewise, antiscarring in between the scarring energies leads to smaller splittings). 
The energy dependence of the distribution of splittings displays EBK quantization, and the 
shape of the smooth envelope is roughly consistent with the prediction of scar theory, though 
the magnitude of the oscillations is not quantitatively predicted by the simple linearized 
dynamics; a better understanding of the shape of the envelope would require extending scar 
theory to the non-linear regime. Also, the distribution of splittings is approximately Porter- 
Thomas with scarring corrections, as we would expect on the basis of scar theory combined 
with the theory of Creagh and Whelan, discussed in Section |. We do not find, however, the 
expected dependence on the Lyapunov exponent of the horizontal periodic orbit. This is 
presumably due to the fact that our calculations do not probe very far into the semiclassical 
limit, our well being only a few wavelengths across in the transverse (y) direction. Finite- ft 
effects cut off the far tail of the splitting distribution at all energies. 

According to Eqs. (0) and @, suitably generalized to the chaotic double- well poten- 
tial in two dimensions as discussed in Section |, the rescaled resonance widths in a single 
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metastable well, the potential of which agreed with the double-well potential we are us- 
ing for x < +xo/y/2, would have the same distribution as the rescaled splittings we have 
computed. Thus, our results imply a non-statistical distribution of resonance widths in a 
chaotic metastable well. In view of its importance for chemical physics, this conclusion 
deserves further investigation. 

Finally, we discuss the prospects for many- dimensional systems. If there exists an un- 
stable periodic orbit near the real continuation of the optimal tunneling path in a double 
well or metastable well, scarring and antiscarring will again play a role. The only question 
is whether the degree of instability is small enough for scarring to be important; for a Lya- 
punov exponent A large compared to unity the short-time envelope approaches the uniform 
limit of random matrix theory. However, as long as the sum of all instability exponents in 
directions transverse to the reaction coordinate does not become large, scarring effects are 
expected to appear, just as in the two-dimensional case discussed in the present paper. 
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FIG. 1. Typical eigenfunctions for the double- well potential in the (a) near-integrable case 
without Gaussian perturbations, E = —11.055, and (b) chaotic case with Gaussian perturbations, 
E = —12.063. Only one side of the well is shown in each case, with the x-axis running horizontally 
from —6 to and the y-axis vertically from —2 to 2. The barrier is located on the right side at 
x = 0. 
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FIG. 2. Contour plot of the potential for a typical member of the ensemble. The perturbation 
at x = ±2, y = ±1 leads to a symmetrical distortion of the contours, which would be more rounded 
in the absence of a perturbation. The contours range from V = —300 near the bottom of the 
potential to V = +200 above the top of the barrier. 
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FIG. 3. Level splitting versus energy E for the 15195 eigenstates between E = and E = —50 
in the ensemble of 625 double- well potentials described in the text. 
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FIG. 4. Rescaled level splitting versus energy E as in Fig. ||. 
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FIG. 5. Rescaled level splitting versus action/27r with eigenstates as in Fig. ||| 
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FIG. 6. Rescaled level splitting versus the overlap of the eigenstate with a Gaussian on 
horizontal periodic orbit, with eigenstates as in Fig. |3[ 
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FIG. 7. Overlap versus energy E, with eigenstates as in Fig. |3]. 
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FIG. 8. Mean rescaled splitting s/{s) versus action/2-7r, data points; short-time envelope from 
scar theory using Eq. @, solid line; short-time envelope from scar theory using Eq. (|5|), dashed 
line. 
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FIG. 9. Distribution of rescaled splittings for numerical data with a = 0.5, solid histogram; 
numerical data with a = 0.63, dashed histogram; Porter-Thomas distribution (without correction 
for finite h, see Fig. Ill]), dashed line. 
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FIG. 11. Distribution of rescaled splittings for the numerical data with a = 0.5, histogram; 
Porter-Thomas distribution with finite-^, correction for N = 6, dashed line. 



22 



